Universal time-evolution of a Rydberg lattice gas with perfect blockade 
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We investigate the dynamics of a strongly interacting spin system that is motivated by current 
experimental reahzations of strongly interacting Rydberg gases in lattices. In particular we are 
interested in the temporal evolution of quantities such as the density of Rydberg atoms and density- 
density correlations when the system is initialized in a fully polarized state without Rydberg ex- 
citations. We show that in the thermodynamic limit the expectation values of these observables 
converge at least logarithmically to universal functions and outline a method to obtain these func- 
tions. We prove that a finite one-dimensional system follows this universal behavior up to a given 
time. The length of this universal time period depends on the actual system size. This shows that 
already the study of small systems allows to make precise predictions about the thermodynamic 
limit provided that the observation time is sufficiently short. We discuss this for various observables 
and for systems with different dimensions, interaction ranges and boundary conditions. 
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I. INTRODUCTION 

Current experimental and theoretical progress in the 
control of ultra cold gases of Rydberg atoms has re- 
vealed remarkable insights into the dynamic and static 
properties of strongly interacting many-body systems 
[Ij [2] . The strong interaction between atoms in highly 
excited (Rydberg) states can be ten orders of magnitude 
larger than the interaction between ground state atoms 
and thus make ensembles of Rydberg atoms particularly 
suited for the study of strongly correlated phenomena. 
The emerging strong correlations are most strikingly re- 
flected in the so-called Rydberg blockade |31 0] - a pro- 
nounced suppression of the probability to excite a Ryd- 
berg atom in the vicinity of another. A number of exper- 
iments have succeeded in demonstrating this excitation 
blockade and provided evidence for the coherent nature of 
the excitation dynamics in an impressive fashion [5l410j. 
From the theoretical side, ensembles of Rydberg atoms 
can often be accurately modeled by spin Hamiltonians. 
This approach has proven to be an extremely valuable 
tool and has delivered insights into the correlated ground 
state phases nTI-[T8] as well as into the dynamics of these 
systems 19 26]. 

A particularly simple spin-model describing the quan- 
tum dynamics of an arrangement of Rydberg atoms was 
explored in Refs. [251 EH EH] by some of the authors. 
There we considered atoms placed on a one-dimensional 
regular lattice - a situation which since very recently can 
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FIG. 1. (Color online) a: General behavior of the time evolu- 
tion of observables such as the density of excited atoms pit) or 
the density-density correlation functions gi(d = 2,t') between 
atoms separated by two sites (lattice with L = 25 sites), b: 
Short time behavior of the density of excited atoms. In the 
upper panel, for lattice sizes L = 10 and 12. In the lower one 
same quantity for L = 16 and 18. One can observe how the 
universal regime grows with L: it is given by t < 3 for L > 10, 
and by i < 7 for L > 16. 
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also be achieved experimentally 1551130]. We moreover as- 
sumed a perfect excitation blockade for neighboring sites, 
i.e., no two neighboring atoms can be simultaneously ex- 
cited to the Rydberg state. For this simplified situa- 



tion, we studied in particular the time-evolution of local 
observables such as the density of Rydberg atoms and 
density-density correlation functions. We found that the 
temporal behavior of all of these observables is qualita- 
tively similar: The short time dynamics is characterized 
by strong oscillations while for long times a steady state 
is established in which all quantities fluctuate with small 
amplitude around a constant value (see Fig. [Tk) . The fo- 
cus of these previous works was on the long-time behav- 
ior, i.e., on the investigation of the origin of this (seem- 
ingly thermal) steady state. Indeed, in Refs. |53J [3T] 
we showed that this steady state can be understood as a 
maximum entropy state and that the evolution into it is 
governed by a Fokker-Planck equation. 

The aim of this work is to study the short-time be- 
havior of observables in a Rydberg lattice gas. In par- 
ticular, we will answer the question of how well the nu- 
merical results obtained for finite systems approximate 
the corresponding results in the thermodynamic limit. 
We find that for sufficiently short times expectation val- 
ues of physical quantities are independent of the system 
size and, hence, universal. An example for that is de- 
picted in Fig. [ijs with the evolution of the density of 
Rydberg atoms. The universal time interval grows with 
increasing system size. We show that there are universal 
functions for the expectation values of general observ- 
ables in the thermodynamic limit to which the results 
of finite systems converge at least logarithmically. This 
means that even the study of finite systems gives access 
to properties of the thermodynamic limit provided that 
the propagation time is sufficiently short. We discuss 
which observables exhibit this convergence to the univer- 
sal behavior, and the validity of these results for systems 
with different geometries, interaction ranges and dimen- 
sions. Note that, based on the so-called Lieb-Robinson 
bounds |32] , the existence of a thermodynamic limit was 
mathematically proven for quantum lattice systems with 
polynomially decaying interactions [33ll35] . In this work, 
we calculate the rate of convergence to this thermody- 
namic limit and show that it is, at least, logarithmic. 
We provide a systematic procedure to obtain these uni- 
versal functions using the particularly symmetric case of 
a ring lattice. 

The paper is structured as follows. In Section |llj we 
describe in detail the system under study and present a 
derivation of the Hamiltonian that governs its dynamics. 



In Section III we discuss the properties of the Hamilto- 
nian and their implications in the time-evolution of the 
expectation value of any general observable. In addition, 
we introduce some of the mathematical tools we will use 
throughout the remainder of the paper. In Section IV we 



focus on the case of a one-dimensional lattice with peri- 
odic boundary conditions and any blockade radius. Here, 
we prove that the time-evolution of the expectation value 
of a set of operators (including the Rydberg density) in 
a finite size system converges at least logarithmically for 
all times to a universal function (i.e., independent of the 
size of the system) in the thermodynamic limit. The tech- 



nique to derive this universal function is also described 
here. In Section Ivj we analyze the same system but 
with open boundary conditions and show that the results 
found in the periodic case are equally valid. Finally, we 
discuss the generalization of these results to higher di- 
mensions in Section \VT\ The conclusions and open prob- 
lems are given in Section [VII| In Appendices \K\ [B| and [C| 
we provide the bulk of the mathematical proof of the ex- 
istence of the universal functions and the convergence to 
those of the finite-size expectation values. In Appendix 
[D| we present the method to create the matrix represen- 
tation of some operators on the one-dimensional linear 
lattice. 



II. THE SYSTEM AND ITS HAMILTONIAN 

The system we are considering is a one-dimensional 
lattice with L sites and spacing a. Each site shall con- 
tain a single atom which occupies the (local) vibrational 
ground state. The ground state \g) and a highly excited 
(Rydberg) state \r) of each atom (considered as two-level 
systems) are coupled resonantly by means of a laser with 
Rabi frequency f2. The atoms in the Rydberg state inter- 
act via the van-der-Waals potential Cg/ |r| , where Cg is 
the van-der-Waals coefficient [36] and |r| is the distance 
between them. Having only two possible states reduces 
the problem to a spin-1/2 system where one can identify 
Iff) — r\r) = \D and \r) = r''' \g) = |t). The annihilation 
and creation operators r^, rl of a Rydberg state at site 



k satisfy thus a spin-1/2 algebra 



[rj,rk\ 



[rhrl] 



rl=rl =0, 






With these operators the Hamiltonian of the system 
reads 



n = nY^ir. + rl) + vY^^, 



k=l 



k^3 



J. 



where nj. = r\.rk counts the number of Rydberg atoms at 
site k and V = Ce/{2a^) represents the nearest neighbor 
interaction strength divided by two to compensate the 
double counting in the sum. 

In the following we will perform two approximations 
which lead us to the perfect blockade model. First, 
we do not consider the full |r|~ -tail of the van-der- 
Waals interaction but rather a cut-off at a radius Rb = 
Af,a, i.e., only atoms on sites contained in the subset 
Afc = {j '■ \k — j\ < Afc and j ^ k} interact with an ex- 
cited atom located at the fc-th site (see Fig. I2|. The 
Hamiltonian can then be rewritten as "H = Hq + Hi^t 
with 



k=l 



+ rl) 






k=ijeAk 



VkjUkTlj, 
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FIG. 2. (Color online) Perfect blockade model. We consider a 
cutoff of the interaction such that only the atoms inside a do- 
main Afc interact with an atom excited on the fc-th site. Two 
different examples of domains are shown: nearest neighbors 
(Afc) and next-nearest neighbors {A'^.) interaction. Moreover, 
we consider the interaction between excited states so strong 
that within the corresponding domain only one atom can be 
excited to the Rydberg state. 



where Vkj = V/\k — j\^. In addition, we assume that 
the interaction energy between atoms inside the radius 
i?h - often referred to as the blockade radius - is much 
larger than the Rabi frequency of the laser, i.e., ft <C 
V/X^. This means in turn that the probability of having 
two atoms simultaneously excited within the set A^ is 
strongly suppressed. These two assumptions can only be 
satisfied if 






(1) 



Note that these conditions are more difficult to meet the 
larger the number of sites inside the blockade radius Rj,. 
In the next step, we make the excitation blockade man- 
ifest by transforming the Hamiltonian into the interac- 
tion picture with the unitary operator U = e~*^'"'* [T5] . 
This results in a new Schrodinger equation with effective 
Hamiltonian 



H 



fe=i 



VkjTl 



'{rk 



J^ 



-2itnk J2i 



Vkj rij 



where we have used that Vkj — Vjk- Since fc ^ A^., all the 
terms contained in the sum over j € Afc commute with 
nfc. Now, we have to calculate explicitly 



U.,- 



'irk + rl)e ' 



aknk „iQfc„t 



'rk, 



with ak = '2tJ2jeAk ^k^nj. Note that n^ = n^, i.e., Uk 
and ruk = l — rik are projectors onto the subspaces with a 
Rydberg and a ground state atom on site fc, respectively. 
Thus, one has 

jeAk jeAfc 

Next, we neglect in this expression the terms that os- 
cillate with frequencies Vkj for all j e A^ which under 
our initial assumption arc much larger than the Rabi fre- 
quency. This is the so-called rotating wave or secular 
approximation, and gives rise to the effective Hamilto- 
nian 



L 

H = nY,irk+rl)Mk, 
fc=i 



(2) 



with Mk — riiGA '^j- This Hamiltonian couples states 
differing in a single Rydberg excitation, but the factor 
Mfe prevents the creation of a Rydberg excitation at site 
k whenever there is another already present at some site 
j e Afc. 

The natural initial state of the system is the one with 
all atoms in the ground state, i.e., |0) = {^j. |g)^, and 
A/fc|0) = |0). The Hamiltonian then only couples |0) to 
states which have eigenvalue 1 with respect to all Mk- 
As a consequence, only states with at most one excita- 
tion inside the neighborhoods {Afcjj.^-^ are accessible by 
a time-evolution under H. 

The generalization of the previous arguments to higher 
dimensional lattices leads to a Hamiltonian of the same 
form as Eq. (2) The higher dimensional version of con- 
dition 1(1)1 is, however, more difficult to satisfy. For in- 
stance, in the case of a blockade radius Rb — a, Eq. ( 
reads V ^ il ^ V/2^, while for a 2-dimensional square 
and a 3-dimensional cubic lattices, it becomes the more 
stringent condition V ^ il ^ V/2^. We would like to 
remark that even when the conditions for the validity 
of the effective Hamiltonian |(2) | are not strictly satisfied, 
such a Hamiltonian provides a first approximation to the 
dynamics of a Rydberg lattice gas which can be corrected 
by perturbation techniques. 

In the next sections, we will mathematically prove that 
the dynamics of a system driven by the Hamiltonian |(2)| 
exhibits some universal features, which do not depend on 
the lattice geometry or the interaction range. 



III. TIME-EVOLUTION OF OBSERVABLES: 
GENERAL CONSIDERATIONS 

Our analysis begins by discussing a property of the sys- 
tem that holds regardless of the geometry of the lattice: 
The operator (—1)", with n = X^fc^-fc being the total 
number of atoms in the Rydberg state, satisfies 

{(-l)"^rfc} = {(-l)"^r^} = 

which in turn yields 



{i-iy\H}^o. 



(3) 



This anticommutation relation implies that H only cou- 
ples states whose total number of Rydberg excitations 
have opposite parity (—1)". Although (—1)" is not a 
symmetry of H, the fact that both operators anticom- 
mute has some implications in the corresponding spec- 
trum and dynamics. 

First, if 1^) is an eigenvector of H with eigenvalue 
E, the state (— 1)"|^) is also an eigenvector, but with 
eigenvalue —E. Therefore, the energy spectrum of H 
is symmetric with respect to -E = 0. Moreover, any 
eigenvector |\I>) with eigenvalue E ^ must be orthog- 
onal to (-1)"!*). Thus, (^|(-1)"|*) ^ 0, which yields 
(^g|*g) = (*o|*o> = 1/2 where |*e) and |*o) are the 
projections of |^) onto the subspaces with even and odd 



total number of Rydberg states, respectively. We con- 
clude from this that any eigenvector of H with non zero 
eigenvalue has equal probability of having an odd or even 
number of Rydberg excitations. 

For some particular choices of the lattice geometry and 
the blockade range A^, the system and its Hamiltonian 
H possess additional symmetries. For instance, in the 
case of a ring lattice (one-dimensional lattice with peri- 
odic boundary conditions) the system is invariant under 
translations and reversal of the sites. Any symmetry of 
the system splits the state space into several i/-invariant 
subspaces (the eigenspaces of the symmetry) , decoupling 
the evolution into lower dimensional ones with the same 
spectral properties described previously for H. If this 
symmetry also commutes with (—1)", the corresponding 
subspaces are (— l)"-invariant too. As a consequence, 
any odd-dimensional subspace which is simultaneously 
invariant for H and (—1)" must enclose a zero energy 
eigenstate of H. Otherwise, the number of orthogonal 
eigenvectors of the reduced Hamiltonian should be even 
due to the symmetry of its spectrum with respect to the 
origin. 

The anticommutation relation |(3)| also bears impli- 
cations for the dynamics. In the Heisenberg picture 
the evolution of any observable A is given by A{t) = 



itH 



At 



itH 



The evolution of this arbitrary observable 



A is an entire function of t, i.e., an analytic function on 
the whole complex plane when t is considered as a com- 
plex variable, so that its power series must converge for 
any value of t. The corresponding power series expansion 
yields 



OO _j 

A{t) = Y^-iH^s.d=^{A), 



(4) 



J=0' 



where we have defined the adjoint endomorphism 

ad^^(A) = [H, A], = [H, [H, A],_i] , ad°,(A) = A. 

Let us assume that [A, (—1)"] = 0, i.e., the operator A 
conserves the parity of the number of excitations. Then, 
the relation [(3)1 yields 

(-l)"^(t)(-l)" = A{-t). 

If |<I>) and 1^) are two eigenstates of n with total number 
of Rydberg states n<s, and n*, respectively, the matrix 
element {^\A{t)\'^) satisfies 



imi-m) = (-1) 



imim)- 



As a consequence, (<I>|^(t)|\l') is an even or an odd func- 
tion of t depending whether n$ and n^ have the same 
parity or not. Setting |<I>) — |^), we find that the ex- 
pectation value of A starting from an initial state with a 
definite number of Rydberg excitations is an even func- 
tion of t. 

We now apply these results to the time-evolution of 
the expectation value of the total number of Rydberg 



atoms n. We consider the initial state |0), i.e., i 
(0|n(i)|0), so that (n(0)) = 0. By using |Eq. (4~ 



n{t)) ^ 
we see 



that only the even terms contribute to the expansion: 



.(i)) = £c,i2.-. 






d'^in)). (5) 



The first term of the sum is proportional to t^ and can 
be calculated using the algebra of r^ and r], as 

ad^(n) - [H,n] - nj^irk - rl)Mk 

k 

ad^(n) = [H, [H, n]] = 2^^ ^(n,- - mk)Mk 

k 

+2"' E E ^"-A - r]r,)M,M,\^ 

k jeAfe 



with M]^\j — Jl;gA i=^j "^'- Taking the expectation value 
in the last expres sion we obtain (ad^(n)) — —2Q^L and 



in the expansion (5) the first coefficient is ci = Q L. 



Therefore, the expectation value of the Rydberg density 
defined as p{t) = {n{t)) / L is an even function of t with 
the universal behavior p{t) — il^t^ + 0{fl'^t^), regardless 
of the dimension or the geometry of the lattice, and for 
any blockade range. 

In the remainder of the paper we will search for such 
a universality (independence of the size of the system) in 
higher terms of the series expansion of pit), as well as for 
expectation values of other operators. To tackle this task, 
we will extensively use the following property: If A and B 
are self-adjoint, a.dg{A) is self-adjoint or anti self-adjoint 
depending on whether j is even or odd. Consequently, to 
obtain a.dg{A) — [B,ad-'Q {A)] we can use the identities 



adi,(A) = <' ^-d-^(A)+h.c. jeven 

B ad^j^\A) - h.c. jodd, 

i.e., the symmetrization or anti-symmetrization of 
Bad-'g (A) for even and odd j, respectively. 

For convenience, we will set from now on the energy 
scale ri = 1. The change t ^ Vlt allows us to recover the 
dependence on Q. in the evolution. 



IV. THE ONE-DIMENSIONAL RING LATTICE 

In this section, we will analyze the dynamics for a ring 
lattice with only nearest neighbor interaction, that is, 
a one-dimensional finite lattice with periodic boundary 
conditions and blockade range K^ = {k — l,k + 1} . This 
system was previously explored in Refs. [Ml HZl [25[|7r] . 

The physical properties of this system will in general 
depend on the number of sites L. We make apparent this 
explicit dependence on L by including a superscript (L). 
The Hamiltonian reads 



H 



[L) _ 



= E^^' 



with 



Hk = TOfc^iAfc ruk+i- (6) 



fe=i 



where the notation A^ — rk + r^, has been introduced. 



According to the discussion in Sec. Ill the Rydberg den- 
sity is given by 



,(^)(t) 



1 



,(i) 



w> = E4'^^^^ 



j=i 



with n'--^' — X]fc=i "-fc! ^^-d 



(2j)! 



where the invariance of the system under translations has 
been used. The periodic boundary conditions imply that 
any index k must be understood mod(L). 

In the previous Section, it has been already shown that 
in general cj^ = 1. We calculate again this coefficient 

using a procedure that allows us to obtain c- beyond 
j = 1. First, we need to compute 

ad^(L)(nfc) = HkUk - h.c, 

where we have used that Hk is the only term of 
ij(^) which does not commute with nj.. This anti- 
symmetrization yields 



ad^(t) (rifc) = mfc_i(rfe - rl)mk+i 



(9) 



In order to compute adjj(L){nk), we first realize that 
among the terms of _ff(^) only -fffc-i, Hk and Hk+i do 
not commute with ad^(i,)(nfe), and hence 



ad^(t)(nfc)= I E Hk+j j mk-i{rk - rl)mi,+i + h.c. 



Upon symmctrization, we obtain 

adject) (rife) - 2a^. + a^. +a^^^;^, 

flfc = mk-2{rk-irl + rl_^rk)mk+i, 
with (5fe = rifc - TOfc. Hence, (ad^(i,)(nfc)) = 2(4) 



and we recover the expected result c-j^ — 1 

We now proceed with the second coefficient Cj 
ad^(i,)(nfc), the same procedure as before yields 



(L) 



(10) 

= -2 
For 



(1) 



,(2) 



.(2) 



adi,,r.,ink)^2A^^>+Ai'>+A^a, 



where 



^i^^ = [E^+.'«i^^]' A^-[T.H..,>^^ 



TABLE I. The most frequent operator 
through the text. 


multiplications ah used 


a\b 


rj 


r] 


Uj rrij 


A. 






mrij 


rrij 
r] r, 



(7) 



Calculating the commutators A^ and A\. and perform- 
(8) ing the anti-symmetrization finally results in 



ad^(., (n,) = 46i.^^ + h'^' + Zh'^l^ + -ih'^' + b^^^^. 



^(4) 






(11) 



&i = mk-i{rk - rl)mk+i 

bk = mk-2mk^iirk - r\)mk+i, 

°k = ^Tik-2{rk-i - r\_^)mkmk+i, 

!.(4) / t t t \ 

As we can see, the level of complexity of the terms 
ad'l,(i,)(nfc) increases very quickly. However, since we are 



only interested in 



„(i) 



it is not necessary to compute 



ad^(i,)(nfe) completely. The only terms of ad^(i)(nfc) 
which contribute to (ad^(i,)(nfc)) are the products of op- 
erators 771/ . Such terms are self-adjoint and, hence, they 
are generated simultaneously with equal coefficients in 
both summands of g^^^ad^ d.) (nfc)-fh.c.. Using the alge- 
bra of rj and A (see Tab. I) the terms of i/'^^)ad^(i,) {uk) 



depending only on m; derive from the products of the 



form Hkh 



(1) 



feOfe 



^(2) 



Hkhl and Hk-ib^. 



(3) 



That is, the only part 



of a,djj(L){nk) contributing to Cj is the symmctrization 



(L) 



of AHkb 



(1) 



k"k 



H.b^^' 



SHk+ib);i, + 3Hk-ib'k' + Hkb^klr 



u(2) 



u(3) 



Hence, we find that {a,dj^(L) (rik)) = —24 and, finally, 



4^^ = ^(ad^(^)K)) 



-1, (12) 

also independent of the size of the lattice L. 



Universality condition of the expansion 
coefficients 



.(i) 



M) 



Neither c\ nor Cj have shown a dependence on the 
size of the ring L. This is, however, only true because 
we have implicitly assumed that L > 2 and thus, the 
boundary effects due to the finite size of the system do 
not play a role. 

Let us explore these finite size effects for the (artificial) 
case of L — 2. The influence of the boun darv bec omes 
apparent in the quantity ad^(2) (ni). While Eq. (9) holds 
identifying mo = m2, i.e., 



i=-i 



i=-2 



ad 



H(2) 



(m) = [Hi,ni] = TOo(ri - r{)TO2 = {ri - r\)m2, 



the calculation of ad^(2) (rii) gives 



TABLE II. For blockade radius Rt = a (A(, = 1), first five 



ad?,(2) (ni) = [iJi + H2,ma{ri - r^mz] (13) universal coefiicients, c,, of the Rydberg density p{t), and of 

pd{t) with d = 2 and 3, Cd^j. For a linear lattice, the first four 
universal coefficients qj. For the blockade radii \b = 2 and 3, 
the first universal coefficients of the Rydberg density pit) are 
also displayed. 



and we see that Eq. (10) is not valid due to two reasons: 
first, some contributions are missing if compared with the 
general expression ad^(L) [ni] = [H^ + Hi + H2, mo{ri — 
r[)m2]; second, in the general case the commutators are 
calculated under the assumption that operators with dif- 
ferent indices commute, but this only holds as long as 
there is no term with indices which are equal mod(L), 
and this is violated by r7T,o(ri — r[)m2 ii L = 2. If we 
proc eed without simplifying mg and m2, the calculation 
(13) oi_adj,(2) (rii) yields only part of the terms obtained 
in|Eq. (10) for the general case. 



(2) 



To calculate explicitly the coefficients c^ ' for L = 2, 

the corresponding ad^(2)(ni) are needed. For j = 1, the 
result is 









Xb-- 


= 1 






\b 


= 2 


\b 


= 3 


j 


Cj 


C2,j 




C3,j 


* 




Cj 




Cj 


1 


1 














1 




1 


2 


-1 


1 




1 


2 
3 




5 
3 




7 
3 


3 


3 

5 


3 
2 




-2 


38 
27 




77 
45 




152 

45 


4 


81 
280 


283 
240 




61 
30 


518 
243 


- 


713 
504 






5 


3023 
25200 


739 
1120 


- 


2393 
1680 


76016 
27207 













ad 



ff(2) 



(ni) = 2Sim2 + rirl + r\r2, 



so that 



„(2) 



= -^(ad?f(2,(ni))=--(-2) = l, 



2! 



2! 



in agreement with the general arguments of |Sec. Ill] 
However, for the next coefficient, we arrive at 



ad' 
ad' 



_H-(2) 



(ni) = 5(ri - r\)m2 + 3TOi(r2 - r\) 



2^ 



_ff(2) 



(rii) — \Q5im2 + Qfni52 + A{rir\ + r{r2). 



.t. 



that yields a different result from Eq. (12) namely, 

2 



.(2) 



{a.d'u^2){ni)) 



1 

4!' 



-16) 



For L = 3, Eqs. |(9)| and (10) remain valid because they 
involve commutators of lT^>= Hi + H2 + H-^ and terms 



with no more than 3 consecutive indices, while Eq. (11) 



„(3) 



~W 



fails. This only ens ures that c\ =1. However, Cj 
— 1, coinciding with Eq. (12) because the perturbation 
created by the periodic boundary conditions does not 
affect the terms of ad^(3)(nfc) contributing to Cj . 

The previous results seem to indicate that ad'l,;^^) (rife) 
and C; are universal for j < L — 1 because they can be 

computed using only the universal algebra of r/ and rj 
without being concerned with the boundary conditions. 
This is indeed true, and the proof is provided in Appendix 
A We denote by Cj the universal value of c, ' for j < 
IT— 1. For illustration purposes, the first five of these 
universal coefficients Cj are given in Tab. II 



B. The universal Rydberg density p{t) 

At least formally, we can define now the power series 

00 



with Cj being the universal expansion coefficients. We 

would naively expect that p^^\t) ^ p{t), but this 

can not be simply implied as a consequence of Cj — 



JL) 



00 c] '■ In fact, for the moment we cannot even 



linii^ 

be sure that the power series in |Eq. (14) has a non null 

radius of convergence. Let us show here that Eq. (14) 

defines an entire function such that p'^)(i) ^ p{t) 

for any value of t. This implies that, although p(t) is 
defined by a power series around the origin, it encodes 
the universality of p'^^^t) in any regime. We follow here 
a complementary approach to that presented in Ref. ^33j 
that will also allow us to calculate the rate of convergence 
to p{t), which will be referred to as the universal Rydberg 
density in the remainder of the paper. 

An analytic expression for all the coefficients Cj and 
c- is not available. However, to prove the above state- 
ments it is enough to obtain bounds |cj|,|Cj- | < hj 
for a set of coefficients bj such that the power series 
X] 7=1 ^j^'' ^^s an infinite radius of convergence. Since 

Er=il^.ll*l''' < TT^ihW: then Y.7=i<^ji^' is abso- 
lutely convergent for any t. Besides, bearing in mind 

that c, = Cj for j < L — 1, we obtain 



00 



W) 



M 



it)\<J2\c,-^'^^\\t\'^<2Y,b,\t\'\ (15) 






which shows that p^^\t) 

Appei 

dents 



— > p(t) at any time t. In 
Appendix iBl we prove that \cj\, \c- \ < bj for the coeffi- 



b- -2-6^^-^-'^^ 
'- (2j)!- 



(16) 



where 



te[o,oo) 



t^- 



In Appendix [C] we show that 



T-j-e 



' hi 7 



and, hence, the radius of convergence of X^i^i ^j^'' i^ 



hm -^ — 

i-i-oo bj 



4lim(2j)(2j-l)^^^ 

- Hm (Inj) = oo. 
9 i— >-oo 



As we pointed out previously, this proves at the same 
time that p{t) is defined for any value of t and that p{t) 
gives at any time t the asymptotic of p^^\t) for large L. 

Convergence rate to the universal behavior. 

Our procedure has the advantage of making it possible 

to estimate the rate of convergence of p*^^-' (i) ^ p{t) . 

From Eqs. |(15)| and |(16)[ we obtain 

\p{t)-p^'\t)\<iy^^6'^t\'^. 



^f-.m 



Appendix [C| proves that 



< 



1 

UJ1UJ2 ■ ■ ■ ^j 



where ujj is the solution of loj + Inwj = 1 + Inj, which 
increases logarithmically with j so that Uj ^ In j. There- 
fore, 



\p{t)-p<^^Ht)\< 



\^ 



62j'|i|2j' 



3 ^ UJ1UJ2 ■ ■ ■ UJ2j 



which means that an upper bound for the error \p{t) — 
p^^\t)\ is given by 



-2LU|2L 



Ei^^it)^l' \'\ (f^ 



e^j'lipj' 



3 Wi • • • (jJ2L 






^ ^2L + 1 ■ ■ ■ '-^2L+2j 



Bearing in mind that uij grows with j, we find that the 
rate of consecutive error bounds satisfies 



E(^\t) 36|t|2 



mt\ 



E(^-^){t) UJ2L-1UJ2L i^oo ln(2L - 1) ln(2L) ' 

i.e, the convergence p(-^^(t) ^ p{t) is, at least, loga- 
rithmic. 

The previous mathematical results suggest the follow- 
ing physical interpretation: For large i, the deviation 
from the universal behavior defined by p{t) is related to 
the finiteness of the lattice. At the beginning of the evo- 
lution, the boundary conditions have little influence on 
the dynamics because the predominant interaction is the 
one between nearby atoms. Mathematically, the first co- 
efficients of the power expansion of {nk{t)) depend only 



on terms Hi with I close to fc, so that the evolution is 
'ignorant' of the boundaries. This is the origin of the 
universality observed for short times in Fig. [TJd. When 
t increases, higher degree terms of (nfc(t)), which include 
contributions from Hi with / far from fc, become signifi- 
cant. This means that the (indirect or mediated) inter- 
actions between distant atoms become more important 
as time passes. Finally, when t is so large that the medi- 
ated interaction between atoms on sites at a distance of 
L or more sites becomes relevant, the universality is lost. 
All this is a consequence of the fact that there is a finite 
speed at which information propagates, which is expected 
from the local interactions that govern the dynamics of 
the Rydberg gas [32] . 

This also means that the actual length of the universal 
time interval shown in Fig. [lb depends on L: The lattice 
size fixes the number of coefficients c, which are uni- 
versal and hence the range of values of t for which p^^' (i) 
will not differ from the universal Rydberg density p(t). 
Indeed, like the number of universal coefficients Cj, the 
range of the universal regime must increase with the size 



(i) for L' >L. 



At this 
^p{t) 



of L, and during it p^^'^{t) « p^^'' 

point, the fact that the convergence p^^\t) 
takes place for any t means that we can just increase L in 
order to extend the length of the universal time interval. 
For finite values of L, these ideal power expansions can 
be calculated numerically in an alternative way using ma- 
trix representations of H^^' and rv-^' . Such computations 
can be simplified by reducing the operators to the maxi- 
mally symmetric subspace under translations and reflec- 
tions (see Ref. [27 ). As an example, we have performed 
this calculation for L = 18 obtaining all the universal 
coefficients until j = 17. A comparison between this re- 
sult and the numerical time-evolution of a system of the 
same size with closed and open boundary conditions is 
presented in |Fig. 3[ Here the universal regime still lies 
in the transient phase and does not extend to the steady 
state regime. This universal time will increase with in- 
creasing number of sites on the lattice and, consequently, 
number of universal coefficients. 




1 1.5 2 

t (units of n"'') 

FIG. 3. (Color online) Comparison between the short-time 
evolution of the Rydberg density for a ring (blue solid line) 
and a linear (red dotted line) lattice with L — 18, and the 
universal Rydberg density X] =1 '^j*^"' (magenta thick-dotted 
line). 



Beyond the nearest neighbors blockade. We will 
consider now a blockade radius Rb larger than the lattice 
spacing a. Suppose that R}, — AfcO, such that the block- 
ade range is Ak = {k — Xf,, . . . ,k — l,k + 1, . . . ,k + Ah} 
for any site k of the ring. This yields an effective Hamil- 
tonian i/*^^^ — J2k=i ^k with 

Hk = TOfc-Ab • • • "T-fe-iAfc mfc+i • • • nik+x,. (17) 



correlation function converges to 1 as t — >■ 0. From 
the universal behavior (nfc(t)) ^ t^, we conclude that 
{nk{t)'ni,^^{t)) ^ t^ should be universal too, so that 

•^d 1 ~ *-* ^^'^ '^d 2 ^ ^ which is the first non-zero con- 



For the Rydberg density, we have shown in Appendix 
Althat ad^(^) (uk) and Cj are independent of L for j < 





TX- 1)/A6. Thus, there exists a universal set of expansion to the case of the Rydberg density ad^^,^, (n^nfe+d) and 



tribution in [Eg. (18) being universal. Indeed, this can 
be shown foll owing sim ilar commutator computations to 
those done in lScc. Illl for the P universal behavior of the 
Rydberg density. 

Concerning the universality of the rest of the coeffi- 
cients c'j , it is proven in Appendix A that, analogously 



coefficients Cj = limL_>oo 4 and a universal Rydberg c^ij are independent of L for j < L-d. This gives rise to 

density p{t) = J2'jLi '^jt'^'' Just hke in the case Aj, = 1. 
For Aft = 2 and A^ = 3, some of the universal coefficients 
Cj are provided in 



a set of universal coefficients Cdj , some of which are given 
Tab. Ill for d — 2 and 3, that allow us to construct the 



Tab. II 



The bounds given in Eq. (16) 
for the coefficients of the Rydberg density, become now 
(see Appendix [b]) , 



^j\' 



|cf^|<2(6Ah)^-'--i^^?^±i^^^ 



im 



which proves that the Rydberg density p*^^^ (t) converges 
for any value of t to the universal Rydberg density p{t) , 
no matter the size of the blockade radius. 

Apart form the first term i^, which is completely 
universal, the universal Rydberg density for a one- 
dimensional lattice only depends on the blockade radius. 



C. Expectation values of other operators 

So far, we have only studied the Rydberg density. How- 
ever, these results can actually be generalized to other lo- 
cal observables. Let us consider, for instance, the density- 
density correlation function 



92{d,t) 



{nk{t)nk+d{t)) 



{nk{t)){7ik+d{t)) 
and, in particular, the expectation value 

pf\t) - {nk{t)nu+d{t)). 

with d > 2, since (rikUk+i) = due to the blockade 
condition. 

As in the case of p'^^^ (t) , p^ ' {t) is an entire function 
even in t and vanishing at i = with power expansion 



JL) 






(18) 



where 






^d,j 



(2.?)! 



(ad^(t)(»^fc"-fc+d))- 



For very short times, atoms on separated sites evolve 
approximately independently. Thus, the density-density 



universal function. 

Some of the previous findings can be extended to any 
observable. To formalize this statement we consider ob- 
servables that are independent of the lattice size L and 
that are constructed as a sum of words. A word is de- 
fined as an operator A = Xk-^x^^ ■ • • a^fc with fci < fc2 < 
••• < kj, built using the letters xi S {r/,rj , m;,?!;}. 
Each of these words is characterized by its length i{A) = 
kj ^ ki + 1 and its number of single letters xi € {n, rj}, 
s{A). Let us consider the expectation value of a word 
A with length £ — i{A) and number of single letters 
s^s{A), 

p^^\t)^{A{t)). 

Since the relation (— 1)"A(— 1)" = {—XyA is accom- 
plished, then (-l)"A(i)(-l)" = (-l)M(-t) at any 
time. Hence, p\^ {—t) — {—iypA{t), i.e., the parity of 
Pj^ (i) is determined by the number of single letters of 

A. As a consequence, the entire function pj^ (t) can be 
written as 



p'A'\t)- 



oo 

i=0 



-A.j'- 






M^WM)) 



where only powers with the same parity as s contribute. 
In this general situation, we can ensure that ad^„(j^) [A) 

and 0^4 , are universal at least for j < {L — £)/2 (see Ap- 



pendix lAl), and we denote by ca^j these universal coeffi- 
cients. Furthermore, the bounds 



|cAj|,| 



'A, J I 



< 12^ 



S+^/2-l 



J 



obtained in Appendix [Bj t ogether with the asymptotics 
of Hj given in Appendix |C| allow us to prove that 

oo 
PA{t) ^^CA^jt^ 

defines an entire function which picks up the universality 
of PA\t)- PA{t) ^"^°°> p^^(i) for any t. 



9 



This result shows that the universahty observed in 
the time evohition of the Rydberg density is not caused 
by the specific form of the corresponding operator, but 
rather constitutes a general feature of this quantum sys- 
tem. 

For a blockade radius larger than the lattice spacing, 
it has been shown in Appendix A that ad'L(j^) {A) and 

0^4 j are universal at least for j < (L — £)/(2Afc), giving 

sense to the universal coefficients c)j '- ^ c^j and the 

universal expectation value PA (i) — J27Li'^A,jtL Finally, 
the bounds 



\CA 



■I Ic^"^^ 



<(12A,)^''MMzi 



uj, 



(19) 



obtained in Appendix [BI show that p]^ (t) ^ pA (t) 

for any t. 

Convergence rate to the universal behavior. 

We have seen that, in a ring under nearest neighbor 
blockade, the convergence of the expectation value of 
the Rydberg density to its universal counterpart is at 
least logarithmic. This estimate is generalized here to 
any other translation invariant or L-independent opera- 
tor A, under any blockade range Af,. Since c)j ' ~ caj for 
j < {L — (.)/{2\), using the bounds [(19)| in this general 
case, we find that 

\pAit)-pf{t)\<2 J2 b,\ty =: E^^Ht). 

In order to obtain the rate of consecutive error bounds 
E^^') and £;(^+2^6) we need to calculate 



'^^' = 12A,. ' 



Kj+e/i2\b) 



b] ""j + lKj+^/(2A,)-l 



From Appendix [C] we get the following inequality 
^<12A,-^'^^ ^ 



] + l UJ 



J+2 



with Wfc being the solution of w^ -f Inw^ = 1 -f In A;. Since 
LOk is increasing in fc, then 



Oj + l 



< 12Af,max < 1, 



L 



L 



2Xh 



1 



J > 



L-l 



2A,, 



which leads us to 



^(L+2A,)(^)^2 



E 



h+iw^' 



< 12Abmax { 1, 



L 



L-e + 2Xh 



-E^^\t). 



Finally, bearing in mind that lo^ ^ Ink, we obtain the 
upper bound for the ratio 



EW(t) 



< l2Xh max < 1 



L-i + 2Xb\ oj 



\t\ 



r^ 12Xb- , 



In 



2At 



which again shows that the convergence is at least loga- 
rithmic. 

In summary, the universality of the evolution of the 
Rydberg density under nearest neighbors blockade also 
holds for the evolution of the expectation value of any 
other operator with an L-independent shape as well as for 
any operator that is invariant under translations, and for 
any blockade radius. In addition, let us emphasize that 
the ring lattice is not merely a useful model to simplify 
the computations due to the symmetries, but also allows 
us to obtain the exact coefficients of the power expansion. 
That is, it grants us access to the universal behavior of 
the system independently of the lattice geometry, as it is 
discussed in the next Section. 



V. LINEAR LATTICE - BOUNDARY EFFECTS 

We will now investigate a one-dimensional lattice with 
open boundary conditions. Considering again nearest 
neighbors blockade, the effective Hamiltonian i/^^-' for 
a lattice with L sites reads as in Eq. (6) but with 
niQ = m^+i = 1. The analysis of the short time dy- 



namics is similar to the ring lattice but the boundary 
effects are more pronounced. The reason is that, con- 
trary to the ring, the atoms located at the edges of the 
line 'notice' the boundary conditions at any time. How- 
ever, the theoretical approach followed previously in case 
of the ring can be extended also to this case. 

We consider the expectation value of the Rydberg den- 
sity p^^ \t) = {rS'"\t) )/L. It c an be expanded as in ex- 



pression 



(7) 



but the 



Eq. (8) 



for the coefficient c- ' of 



t^^ is no longer valid because the linear lattice has no 
translation invariance. Only the invariance under reflec- 
tion remains, so that the contributions of the sites k and 
L — k coincide. Therefore, in the coefficients 



SL) 



1 (-IP 



L (2j)! 



iV ^ 



(ad|(^)(nA;)), 



fe=i 



the contribution of each site k must be computed inde- 
pendently. 

The first coefficient c\ is proportional to the expec- 
tation value (ad^(i,)(nfc)). First, we need ad^(i,)(nfc) — 
[7Jfc,nfe], and to calculate it we follow the same route as 
in the ring when k ^ 1, L, since for k = 1 and L, one has 
ad^(i,)(ni) = Aim2 and ad^(t)(ni) = m^^iA^, respec- 
tively. Generalizing this for any k, it yields 



ad}^(i,)(nfc) = mk-iAkruk^ 



mo = jriL+i = 1. 
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In the same spirit, one can write the operator 

ad^(^,(nfc) = [g(fc-i) + g(fc) + Jj-(fc+i),ad^ff„( nfc)] with 



"^o = ^TT-L+i — Ij and rn_i = 771^+2 = 0, so that Eq. (10) 
remains valid for any k. From here we obtain again that 
(adjLf(t) {nk)) = 2((5/c) = —2 for any k and 



JL) 



^X:(ad^(.,K)) = ^(-2L) = l. 



2L 



fc=i 



2L 



This result was expected because we know that the i^ 
behavior of the expectation value for the Rydberg density 
is independent of the geometry of the lattice. 

The previous arguments show that the coefficients of 
the power series expansion of the Rydberg density on the 
linear lattice can be inferred from the ad-'„ fj^^{nk) on the 
ring by setting rriQ = m^+i — 1 and ruk = if fc < or 
k > L + 1. This rule not only simplifies the calculation 
but at the same time provides a direct comparison with 
the ring and its universal Rydberg density. Following this 
recipe we find that the next two coefficients are 



„W 



1 



2 



and 



r,W 



1 



38 
27Z 



Hence, except for c\ , the coefficients c, are not uni- 
versal, even for large L. However, for L sufficiently large, 
all sites give the same contribution to the j-th coeffi- 
cient, except for the j — 1 sites counted from either end 
of the lattice. For each of these boundary sites, a fraction 
Pi, . ■ . ,Pj-i of the universal contribution is missing due 
to the condition ruk = ii k < or k > L + 1. Hence, 
we can rewrite 



„(i) 



i^' 



2pi 



2P,-1 



1 
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L 



where Qj — 2(pi + • • • + pj_i). Although the coefficients 

Cj of the open lattice depend on its size L, the values 
of Qj are i-independent provided that L is large enough 
compared with j, in particular if j < L + 1. This re- 
sult can be checked with an alternative computation of 
ad^(j^) (n^^') for concrete values of L, using the matrix 
representation of i/*^-^^ and n^^-' developed in Appendix 
[D| The first universal values of qj obtained with these 
calculations are given in |Tab. I"I| 

In general, there is no strict universality in the coef- 
ficients of the power expansion for p(^)(i) on the line. 

Nevertheless, as for a ring lattice, c- ^ Cj because 

Cj — Cj + 0{\/L). Furthermore, p^^\t) also converges 
for any t to the universal Rydberg density p{t). The 
reason is that this convergence in the ring is simply a 



consequence of a bound for c, given by an overestima- 

tion of the number of words in ad^,^) (n^) contributing to 

c, (see Appendix bL Since the linear lattice reduces the 
number of such words due to the conditions m^ — for 
A; < and k > L + 1, the bound holds here too. As a con- 
sequence, we also expect the presence of a universal short 



time regime with a range that increases with the size of 
the system. A comparison between the Rydberg density 
for a L = 18 linear lattice and the universal Rydberg 



density including terms with j < 17, i.e., ^ 



17 



jf^, is 



presented in Fig. 3 



These results can be generalized to other operators and 
other blockade ranges. 



VI. MULTIDIMENSIONAL LATTICES 



For completeness, we provide some ideas on the gener- 
alization of these results to systems of higher dimensions. 
The analogue of a one-dimensional ring lattice of size L in 
higher dimensions is a toroidal lattice defined by the peri- 
odic boundary conditions k+Lej = k for any spatial basis 

vector Cj. The coefficients c, of the Rydberg density 
will depend not only on the blockade radius but also on 
the dimension of the system. Higher dimensional lattices 
yield different quantitative results for non-integer values 
of Af, due to the presence of neighbors at non-integer dis- 
tances, absent in the one-dimensional case. Neverthe- 
less, the short time behavior will be qualitatively similar: 
The first Taylor coefficients of the expectation value of 
a (L-independent or translationally invariant) operator 
are universal for L large enough because the correspond- 
ing commutators involve sites with a relative distance 
not greater than L. Furthermore, the number of uni- 
versal coefficients increases linearly with L. As in the 
one-dimensional case, these universality properties are af- 
fected by non-periodic boundary conditions (e.g., square 
and cubic lattices): The coefficients differ from the uni- 
versal ones by 0{1/L) terms because, for large L, the 
number of boundary sites is of order 1/L compared with 
the total number of sites of the lattice. 

Hence, we can at least formally define a power series 
expansion that represents the universal time evolution of 
expectation values of operators such as the Rydberg den- 
sity. However, the reasoning used earlier in this paper to 
prove that this series actually defines a function (i.e., that 
the power series with the universal coefficients is conver- 
gent) does not apply to dimensions larger than one. The 
cornerstone of such methods is the bound obtained for 
the Taylor coefficients by overestimating the words gener- 
ated when applying ad'^;^) to the corresponding operator 
(see Appendix mn . The number of such words grows ex- 
ponentially with the dimension, which renders the tech- 
niques useless in dimension 2 or higher. This, of course, 
does not mean that there exists no such function as the 
universal Rydberg density in higher dimensions. In fact, 
its existence was proven previously in Ref. j33j although 
the method to calculate the corresponding universal func- 
tion remains, to the best of our knowledge, unknown. 
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VII. CONCLUSIONS 



Appendix A: Universality on the ring 



In this work we have shown that expectation values 
of observables in a one-dimensional Rydberg lattice gas 
converge to universal time dependent functions in the 
thermodynamic limit. This convergence takes place for 
any time, for any size-independent or translationally in- 
variant observable, for any blockade radius, and for open 
and periodic boundary conditions. From a mathematical 
point of view, the universality can be understood as an 
asymptotic behavior when the size of the lattice goes to 
infinity. The expectation values of the considered opera- 
tors converge (at least logarithmically) to some universal 
ones which could be associated to a model living in an 
infinite lattice. We have also shown how to obtain the ex- 
pansion coefficients of these universal functions by inves- 
tigating the dynamics on a ring lattice. The difficulty in 
the analysis of this asymptotics stems from the fact that 
such infinite lattice model is not well defined since the 
corresponding Hamiltonian has no meaning as an opera- 
tor on a Hilbert space. Higher dimensional lattices also 
share in general these properties and the expectation val- 
ues actually converge to universal ones as the lattice size 
increases. The universal expectation values given this 
asymptotics depend only on the operator at hand, the 
blockade radius, and the dimension of the lattice, but 
not on its particular geometry. 

From the practical point of view the results presented 
in this paper mean that even the simulation of small sys- 
tems allow to gain accurate insights into the behavior 
of macroscopic many-body systems, provided that the 
simulation time is sufficiently short, so that boundary 
effects are negligible. This has recently been exploited 
in the Ref. [35j which discusses the simulation of open 
quantum systems with local system-bath coupling. 
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1. Rydberg density 



We will prove here that ad'^(^) (uk) and the coefficients 

Cj of the Rydberg density are universal for j < L — 1. 
The proof of these results relies on the following basic 
facts: 

(PI) For 1 < j < L — 1, ad'L(j^)(nfe) is a sum of words 
X = XpXp+i ■ ■ ■ Xp^N which are products of letters xi € 
{ri,rj ,mi,ni}. The initial and final letters are always 
Xp = rUp and x^+at = mp^^. The length l(x) = A^ + 1 of 
each word x is not greater than j -I- 2 and p < k < p + N. 

(P2) We denote by s{x) the number of single letters 
xi € {fhri} in the word x — XpXp^i ■ ■ ■ Xp^^. If 2 < 
i{x) < L and Xp = mp, Xpj^N — "m-p-^-N, the commutator 
[H^-'"\x] is a sum of words x' with £{x') > £{x) and 
s{x') — s{x) ± 1. Furthermore, 



l{x') > £{x) ^ i{x') = i{x) + 1, s{x') = s{x) + 1. 

Proof. We will prove (PI) and (P2) by induction on j. 
Equations 1(9) I and 1(10) [ show that they hold for j = I and 
L > 2, and j = 2 with L > 3, respectively. For j = 1, 
ad^(i,)(nfe) is composed by two words, with £{x) = 3 
and s{x) = 1. For j = 2, we find in ad^(i,)(7ifc) words 
with £{x) — 3 and s{x) — 0, and others with £{x) — 4 
and s{x) — 2. Assuming that they are true for an index 
j — 1 < i — 2, we will prove that they also hold for 
j. Under the induction hypothesis, ad''„^,^^{nk) should 
be a sum of the words obtained from [iJ*^-^^ , x] , where 
X — rripXp+i ■ ■ ■ Xp+N-iTTT-p+N are words of aA''~^L^(n]^) , 
which must satisfy p < k < p + N and £{x) < j -I- 1. 

Since £{x) < j + 1 < L, the word x contains at 
most as many letters as terms Hi are in H^^' . Hence, 



[7I(^), 



[E 



P+N 



Hi , x\ where the sum runs over dif- 



ferent indices mod(L). This fact, together with Tab. I 
yields the decomposition 

Y7i=p Hix = mp^iApmpmp+iXp+iXp+2 ■ ■ ■ Xp+N-irrip+N 
+ rrip f I]r=p+T^ Hi) ^P+i ■ ■ ■ Xp+N-irrip+N 

+ rUpXp+i ■ ■ ■ Xp+N'-2'mp+N-lXp+N-l 

'Ap+N'nip+N'mp+N+i 

= mp^irlxp+iXp+2 ■ ■ ■ Xp+N-irrip+N 



{T.y 



Vp+l ■ ■ ■ Vp+N-l I "Ip+W 



4 



-I- nipXp+i ■ ■ ■ Xpj^N^2Xp+N-irpj^iq'mp+isi+i- 

where xi = ?n;X/ are new letters and ^ is a sum over 
the words y = j/p+i • ■ • yp+N-i- This equation shows that 
^^^ Hix is a. sum of words x' with operators mi at the 
right and left extremes. The words have length £{x') = 
£{x) or £{x') = £{x) + 1, and thus £{x') < j + 2. Since 
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we have chosen j — 1 < L — 2, the n j < L — 1. For 
J2f=p ^Hi, an analogous relation to (Al) satisfying the 
same properties is obtained. Therefore, the statement 
(PI) is proven. 

In addition, the multiplication by the factors mi does 
not alter the number of single letters but the words xi and 
Aix j alway s differ by one in the number of single letters 
(see Tab. I). Hence, J2f=p ^i^ i^ ^ ^^m of words x' with 
s{x') = s{x) ± 1. Furthermore, if i{x') = t{x) + 1, then 
s{x') — s{x) + 1. A similar analysis for YTi=p ^^i gives 
rise to the same results, and therefore (P2) is proven. D 

a. Universality of adi, ^j^^{nk) for j < L — 1. 
Result (PI) ensures that, for j < L, ad-'^f^^j{nk) — 

[H^-'^\a.d■'„^^y{nk)] involves only commutators of H^-'^^ 
with words x = TTipXp+i • • • a;p-)_Ar_i?7ip-|_Ar of length not 
greater than j + 1. Therefore, for j < L — 1, ad^j^, (uk) 
is a sum of terms [X]f=p ^i ' ^] with different indices I 
mod(L), that makes the calculation universal since it is 
determined only by the universal algebra of r; and rj . 

For j > L, if we compute the commutators without 
simplifying the factors a t the extr emes with the same 



each step x 



(9) 



[;(ff+i) with g > 1 can increase the 



length of the word at most by one unit, and, if this hap- 
pens, the number of single letters must increase by one 
unit. In the first step x^°^ — >■ x^^\ the length is increased 
by two units and the number of single letters by one unit. 
Since x^'^' — n^, the history of a word x = x^"^^' which 
is a product of m; operators must have the same number 
of steps rising and lowering the number of single letters, 
given that s{x) = 0. Hence, at most j steps of the history 
of X can increase the length and, as a consequence, the 
intermediate words x^^' must accomplish i{x^3') < J + 2. 
Let us show by induction on g that the history of such a 
word X must be universal for j < L — 1: We assume that 
x^^^ is a universal word of ad^(i)(nfc). Since i{x^^^) < 
j + 1, analogously to the discussion of the universality of 
ad^ (j^) (rife), we conclude that [_ff^^\x*^s)] is universal for 

j < L — 1. As a consequence, the words x of ad^(j^)(nfc) 
contributing to (adr^,^^) (n^;)) are universal for j < L — 1, 

and thus c, becomes universal for j < L — 1. 
J J — 



2. Other operators 



index mod(L), i.e., as in Eq. (13) ad'L(j^)(nfe) becomes 



a sum of words already present m the universal case, The analysis of {nu[t)nu+i{t)) requires to consider 



but with some of them missing. This is due to the fact 
that the words x — Xp- ■ ■ Xp+N from such a computation 
can have letters with indices which are equal mod(L). 
Then, when doing [Hi^x], in Hix, Hi acts from the left 
on the first factors a;;j_ia;ijX/j+i such that li — I mod(L), 
whereas in xHi , Hi acts from the right on the first factors 
xi^-ixi^xi^+i such that I2 = I mod(L). That is, [Hi^x] 
reproduces the terms Hi-^x and —xHi^ of the universal 
computation of [Hi^ , x] and [Hi^ , a;], but the terms —xHi^ 
and Hi^x are missing, together with the full commutators 
\Hi: , x] for any index /' = / mod(L) such that Zi < I' < l2- 

b. Universality of c- for j < L — 1. As a con- 
sequence of (P2), ad^(i)(nfc) is a sum of words whose 
number of single letters has the same parity as j. The 
only words of adi,(j^) (rife) that give a non-zero contribu- 
tion to (ad^(j^) ("•*;)) are products of TOj-operators. Thus, 
for those words s(x) — and, hence, they appear only 
when j is even. 

Any word x^^^' of adji^^^{nk) comes from a previous 

one a;'^^^""'^' of ad^(~^) (nfc) when doing [H^-'^\x'^^^~^''] and, 

analogously, a;^^^"^^ comes from a;'^^^^^-' of ad^^j^y{nk)- 

Following this procedure, each word x'^-'' of adji^j^^{nk) 
has a history, 



nk 



.(0) 



r.(l) 



,(2j-l) 



,(2j) 



where x^^' , at the generation g, is the word of ad?,(j^) (rife) 

giving rise to a:^^"*"^' (among other words) through 
[ij(^),a;(9)]. 

Result (P2) states that i{x'^9+i)) > (.{x'^a)) and 
5(^(3+1)) = s(a;(9)) ± 1. In addition, (P2) ensures that 



words X — Xi^xi^ • ■ ■ Xi-, Xi £ {r;, r J , ttt,/ , n/ } , with indices 
h < I2 < • ■ ■ < k that could be non consecutive. We 
define their length l{x) by £{x) = k — h +1. 

At each step ad-^f^,, (n^nfe+d) ^ ad-^,^,(nfcnfe+d), the 
length of the words increases by one unit or remains 
invariant, and the number of single letters increases or 
decreases by one unit. In particular, the first step al- 
ways increases the length, and, at any step, an enhance- 
ment of the length also implies an increase of the num- 
ber of single letters. Hence, we obtain the following re- 
sults for j < L — d: i) the operator ad^(i) {ukrik+d) is 
a sum of words with initial and final letters of type mi 
or ni, and with length between rf + 1 and j -I- d -I- 1; 
ii) if j ^ 1; the minimum length is d + 2. Hence, the 
same reasoning followed for the Rydberg density proves 
that ad'^(j^) {nkrik+d) and c]j are independent of L for 
j<L-d. 

We analyze now the general case of an operator A 
which is a word of length i. All that we can ensure is 
that in each step ad-'~^j^^{A) — > ad^(j^)(^) preserves the 
length of the words or increases it by one or two units, 
while the number of single letters increases or decreases 
by one unit. From here and following similar arguments 
as in the previous cases, we find that ad-L^j^, (A) and c)^ ■ 
are universal at least for j < [L — ^)/2. 



3. Beyond nearest neighbors blockade 



Now, we consider a blockade radius Ri, = Af,a with 



the Hamiltonian given in Eq. (17) Concerning the Ryd- 
berg density, if j < {L — 1)/Ah and j > 1, the operator 
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ad^„(j^) (rifc) is a sum of words with left and right ends of 
type mi^imi^2 • • ■ i^i+Xb ^-^id with length between 2Ab + l 
and Af,(j + 1) + 1. This is due to the following facts: i) 
each step ad-^i)(nfe) — >• ad^,;^, (rifc) increases the length 
of the words by at most A;, units, except for the first 
one which increases it by 2Af, units; ii) at any step, the 
number of single letters increases or decreases by one 
unit, and an increase of the length also implies an in- 
crease of the number of single letters. From this, simi- 
larly to the case of nearest neighbors blockade (Afc = 1), 
we find that ad^„(j^)(nfc) and c, are independent of L 

iOTJ<iL-l)/\b. 

For an arbitrary operator A of length £, each step 
ad^f,) {A) — > ad-^(i) (A) increases its length by at most 
2Xi, units, while increasing or decreasing the number of 
single letters by one. This implies that ad-L,j^) (A) and 

0^4 '■ are universal for j < {L — £)/{2Xfj). 



Appendix B: Bounds for c- and c)^ ■ 



1. The case of the Rydberg density 

A very rough bound for \c- \ is given by l/{2j)\ times 
the total number of words of ad^^^^j (n^) because each one 
of its terms contributes to (ad^(^) {'^k)) with the value 1, 
-1 or 0. Since the number of words of ad-^,^, (rife) is non 
decreasing with L and constant for j < L— 1, we conclude 
that 



„(i) 



'-JIH 



< 



W2j 

(2j)!' 



with Wj being the number of words in ad^(j^)(nfe) for 
]<L-l. 

In order to count the number of words in ad'L(j^)(nfe), 
we start from ad^(i,) (nfc), that has only 2 words of length 
3, and make j — 1 steps to arrive at ad'l,(^)(nfe). Ac- 



cording to Eq. (Al) for any word x of ad?,(j^) (rife), the 
commutator [H'--^^,x] gives at most 4 words with length 
e{x) + 1 and 2{i{x) - 2) with length i{x). Thus, if we 
consider only the histories with exactly I steps increasing 
the length of the words, there are ("'7^) possibilities for 
choosing them. Given a certain choice of these / steps, 
the number of histories is bounded by 2-4'(2(Z + l))-'^^^': 
Each of the I steps increasing the length gives at most 4 
new words for each old one; the remaining j — 1 — I steps 
must preserve the length and any of them gives at most 
2(^—2) new words for each old one of length £. This result 
is obtained taking into account that / steps increasing the 
length give a word of length 1 + 3^ so at any intermediate 
step ^ < / -I- 3. 



Therefore, defining kj — maxfgjo.oo) t-' *, one obtains 



1=0 

1=0 



J-1 



w,<J2i\ M2.4'(2G + l)r'-i 



^ ^U^2^-^-^ = 2-&~'^K,. 



(Bl) 



Using this result, we finally find that 

which is the bound for the coefficients we were looking 
for. 



2. The case of a general word A 

Let us assume that A is a word of length £. We split 
the j steps of A — ^ a<i'„^^^{A) into histories with li and 
I2 steps increasing the length in one and two units, re- 
spectively. Then, the number of words wa.j of ad^„(j^) (A) 
can be bounded by 

WA,j< J2 (/, )4'^+'H2(/i+2/2+£-2))-'"-'i-'= 
h+h<j ^ ^ '^' 

^ E fAV-''(^i+^2 + i-ir'-'^<i2-''.,+|_i. 



ii^ii<] 



hk 



The multinomial coefficient (j ■'j ) counts the number of 
ways of choosing /fc steps increasing the length in fc units 
among the j steps in A ^- ad-t.(^) i^A). Starting at A with 
length £, each of the Ik steps generates at most 4 new 
words from each old one. The remaining j — h — I2 steps 
keep invariant the length and, thus, none of them yields 
more than 2{£' — 2) new words for each old one of length 
£'. The bound follows from the fact that £' <£ + li + 2l2. 
Since (ad^(j^) (A)) is bounded by the total number of 

words of ad-L,j^)(A), one finally obtains 



\ca.j\,\ca;j\ < 



(i) I ^ w^ < 12J 



i'^j+e/2-i 



J 



as the bounds for the coefficients. 



3. Non nearest neighbors blockade 

For a blockade radius Rt = A^a, the bound for 
the number of words Wj of ad^(^)(nfe) that generalizes 
Eq. (Bl)lyields 



i-i 



w. 



^E(^ /)2(4A,y(2(Ab; + l)) 
< 2{6Xty-^K^+.,,. 



3-1-1 
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Among the j — 1 steps in ad^(i) (uk) -^ ad-^jj^) (n^) there 
are ("'7 ) ways of choosing I steps with increasing length. 
There are 2 words of length 2Ab + 1 in ad^(L) (f^fc)• Each 
of the / steps increasing the length gives at most 4Ab new 
words from each old one. The remaining j — 1 — ? steps do 
not increase the length and, hence, none of them yields 
more than 2{i—2Xb) new words for each old one of length 
£. This leads to the above inequality because i < Xbl + 
(2Ab + 1) when only I steps increase the length. Then, 
the coefficients of the Rydberg density are bounded by 

|c,i|cf)|<2(6A,)^-i^!i?^ggpi. 

For a word A of length £, a bound for the num- 
ber of words WA.j of adi^(i) (A) follows by splitting the 

j steps of j4 — > ad'L(j^) (A) into histories combining Ik 
steps increasing the length of the words in k units with 
k = 1,2,..., 2A6. This yields the following bound for 

WA,], 



Y,lk<3 ^ ' ^ 



4?''(2(Efc^^fe+^-2Ab))' ?'', 



and using the inequality Ylk klk+i—^Xt < 2Ab( X]fc h + 
^/2A(, — l) , we arrive at 



WA.J < {l2XbyK,^i_^ 



J+2xr-i' 



and hence, 



CA,, I, |4^] I <(12At 






Appendix C: Asymptotics of Kj 

Let a be any positive real number and 

K{a)^ max /,(t), f^{t)=e-'. 

te[o,oo) 



We will prove here that, for any e. 



K{a — e) 
K{a) 



a— voo 



In a 
a 



First, fa{t) is a positive and C°° function in (0,oo) 
such that limf_5.o fait) = limt->oo fa(t) = 0. This implies 
that fa(t) reaches its maximum at a point t e (0,oo) 
such that /^(t = t) = fair) {a /t — 1 — In r) =0, i.e., 



with cr{a) ~ r(lnr)^ and 

K{a — e) 1 



K{a) 



g(j(a) — o-(a— e) 



(C2) 



Both T{a) and a{a) are C°° functions of a € (0, oo) 
because ga{t) is C°° for a,i € (0,oo) and g'^it) < in 
this interval. In addition, 

a' (a) = r' In r(2 + In r) = In r, 

are positive functions for a > and 1, respectively, hence 
T{a) and cr(a) are increasing in the corresponding inter- 
vals: T'{a) = T/{a + t) > for a > because r > in 
that interval by construction; as a consequence, if a > 1 
then T{a) > t(1) = 1 and a' (a) > Inr(l) = 0. 
The mean value theorem yields 

a{a) -a{a- e) = scr'i^a) = elnr(^Q), ^a e (a-£,a), 

so that, in |Eq. (C2)[ 

K{a — s) 1 



«(«) {r{^a)y 



(C3) 



Using Eq. (CI) we find that uj{a) = a/T{a) satisfies 
cj + Inw = 1 + Ina, 



-> oo and 



which proves that Lo{a) is increasing, Lo{a) 

that in th at limit w(a) ~ Ci;(a)-|-lna;(a) ^ Ina. Therefore, 

Eq. (C3) I yields 



K{a-e) fuji^a) 



K{a) 



L 



In a 
a 



where we have used that lima_>.oo(a/Ca) = 1- 

In addition, since r^a) is increasing and ^^ < f^i one 
obtains that 



K(a — e) 1 

> 



K{a) iTia)y 



Lu{a) 



so that, finally, given that t(1) — k(1) — 1 and denoting 
Kn — K{n) and aj„ = a;(n) with n = 1,2,3,..., 



n 



l^n l^n — l 1^2 

Kji — • • • < 

^n-l Kn-2 1^1 WiaJ2 ■ ■ ■ <^n 



Appendix D: Matrix representations on the linear 
lattice 



l-lnr^O. 



(CI) 



The solution of this equation is unique for any a because 
gait) = a/t — 1 — In t is decreasing in (0, oo). Then, 



K(a) = fair) 



^Tlnr _ g(T(a) 



Consider a linear lattice of L sites under nearest neigh- 
bors blockade. The Hilbert state space has dimension 2^ 
and a basis is given by the eigenvectors with Rydberg 
excitation numbers rii, . . . , n^,, which we label with the 
corresponding eigenvalues. The absence of translational 
symmetry leaves only the reflection symmetry to reduce 
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the degrees of freedom, but it is unable to eliminate even 
half of them. Instead of this, we will reduce the degrees 
of freedom working in the blockade subspace, spanned by 
those states without consecutive Rydberg excitations. A 
basis i?*^^^ of the blockade subspace can be inductively 
constructed on L as 



SW = {|0),|1)}, i?(2)={|00),|10),|01)}. 



(Dl) 



where and 1 represent an atom in the ground and Ryd- 
berg state, respectively. The dimension d'^' of the block- 
ade subspace satisfies 

which defines the Fibonacci sequence, given explicitly by 



d(^) 



75 



Hence, the restriction to the blockade subspace reduces 



75 



the degrees of freedom from 2^ to d'^^ 

In the basis B^^^, the matrix representation A/"^^-* of 
the total number of Rydberg excitations n*-^' — X]fc=i "-fe 
is given by the recursion 



^(L)^ 



A/-(i) 



V(^- 







AA(^-2) _^ J 



L-2 , 







/o o\ 



,0 1 



, AA(2) 



1 

Vo 1/ 

where J^ stands for the identity matrix of order dS"^' . 



Concerning the Hamiltonian H^^\ its matrix repre- 
sentation T-L^^^ in the basis B*^^) is defined as follows: 
i?^^-' transforms a state of i?^-^^ into a sum of those 
states obtained changing a 1 by 0, or a by 1 when- 
ever there is n o 1 in the contiguous sites. Bearing in 
mind 



Eg. (Dl) 



the action of H^^^ on i?^-^) is determined 
by its action on B(^-i) (g) |0) and B^^^^) (g) |01). Writing 
B(^-i) g) |0) = B(^-2) (g) |00) U B(-^~3) (g |010) shows that 

On the other hand, 

1 00) are the first d^^-^^ 



where we remind that -B'^ ^^ 
states of B(^). 

The previous identities translate into a recurrence for 
■^(i) given by 

with Oj J- being the null d'-*^ x d'^^-' matrix. This generates 
the matrix Ti.^^' for any L starting from 




H« = 




H(2) = 



/o 1 A 

1 

V 0/ 



Since lO) is the first element of the basis B^'"^ the 



itH 






itm 



expectation value p'^^^t) = i;(0|e 

the (l,l)-coefRcient of the matrix j- e**^ VV^^^ ( 
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